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Abstract: Higgs boson production in association with a hard central photon and two 
forward tagging jets is expected to provide valuable information on Higgs boson couplings 
in a range where it is difficult to disentangle weak-boson fusion processes from large QCD 
backgrounds. We present next-to-leading order QCD corrections to Higgs production in 
association with a photon via weak-boson fusion at a hadron collider in the form of a flexible 
parton-level Monte Carlo program. The QCD corrections to integrated cross sections are 
found to be small for experimentally relevant selection cuts, while the shape of kinematic 
distributions can be distorted by up to 20% in some regions of phase space. Residual scale 
uncertainties at next-to-leading order are at the few-percent level. 



Keywords: [NLO Computations, Higgs Physics, Standard Model, QCD| . 



Contents 



|l]. Introduction ^ 

|2|. Details of the calculation § 
|2.1| Tree-level calculation and approximations ^ 



2.2 Calculation of the NLO-QCD corrections 



U Checks I 
|3|. Numerical results ^ 
EL Summary and conclusions 12 



1. Introduction 

The start-up of the CERN Large Hadron Collider (LHC) marks a new era of high energy 
particle physics. A major goal of the LHC is the discovery of the Standard Model (SM) 
Higgs boson and the determination of its properties |]|, In this context, weak-boson fu- 
sion (WBF) processes have been identified as an important class of reactions. In particular, 
Higgs production via WBF, i.e. the electroweak (EW) reaction qq — )• qqH, where the decay 
products of the Higgs boson in the central-rapidity range are detected in association with 
two tagging jets of large invariant mass, provides a promising discovery channel for the 
Higgs boson ^ ^, |6|. Once the Higgs boson has been found and its mass determined, 
WBF will allow for a determination of its CP properties ^] and couplings to gauge 
bosons and fermions 10|. Combining information from the H — )• rr, H — )• W~^W~ , 



H — )• 77, and H — )• invisible channels, the couplings of the Higgs boson to the top quark, 
tau lepton, and the weak gauge bosons can be constrained with an accuracy dictated by 
the amount of data available. Because of challenging requirements on the ATLAS and CMS 
triggers and large QCD backgrounds, however, the determination of the Hbb coupling in 
Higgs production via WBF remains difficult [ pT[] . Therefore, new search strategies have 
been suggested, such as making use of the sub-structure of so-called "fat jets", resulting 
from the bottom quarks into which a boosted Higgs boson decays in WH and ZH [^] or 
in tiH production [^] at the LHC. Additional constraints on the bottom quark Yukawa 
coupling could be provided by a future high-energy lepton-hadron collider such as the 
CERN Large Hadron electron Collider (LHeC), which offers a cleaner environment than a 
hadron-hadron collider [0, |l5|, 16, [l7[| . 



Alternatively, extra gauge boson radiation in WBF can serve as a valuable tool for 
improving the signal-to-background (S/B) ratio of the Hbb mode in a hadron-collider envi- 



ronment. In Ref. 1 18 1, WH production via WBF has been found to allow for distinguishing 
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a Hbb coupling compatible with the SM from certain other scenarios. The event statistics 
of the WBF WH mode is limited, however, by the requirement of a leptonic decay of the 
W boson |18, 19|. This loss in statistics can be avoided by requesting a hard photon rather 
than a massive gauge boson in association with the Higgs boson produced via WBF |2^, |l| . 
Indeed, it has been shown that the presence of a central photon can improve triggering 
efficiencies for the multi-jet final state needed to select pp — )• -ff (— )• bb)jj events. Moreover, 
due to a large gluonic component, QCD backgrounds to the bbjj final state are less active in 
radiating a photon of large transverse momentum than the quark-dominated WBF signal. 
Additional interference effects have been found to suppress backgrounds even further [pO|| . 
The consequence of the central photon requirement is thus a pronounced increase in the 
S/B ratio, making the channel pp — )• H'jjj a process particularly worthwhile to investigate. 

In view of the importance of this channel, precise predictions for the signal process are 
essential. We provide next-to-leading order QCD corrections to H'jjj production via WBF 
at a hadron collider in the form of a parton-level Monte Carlo program, structured similarly 
to existing code for WBF-type reactions |^^. The program allows for the calculation of 
cross sections and distributions within experimentally relevant selection cuts. In order to 
obtain infrared-safe predictions for the Hjjj final state, we employ the photon-isolation 



criterion of Frixione [23| 



We start with a brief description of the calculation in Sec. In Sec. ^ we provide a 
detailed phenomenological study of H'jjj production via WBF at NLO-QCD accuracy. We 
estimate the theoretical uncertainties of our predictions by analyzing the scale dependences 
of integrated cross sections and present kinematic distributions within different experimen- 
tal settings. The impact of the QCD corrections on various observables is quantified. Our 
conclusions are given in Sec. ^ 

2. Details of the calculation 

2.1 Tree-level calculation and approximations 

At proton-proton colliders, Higgs boson production in association with a photon in WBF 
mainly proceeds via quark-quark scattering processes, qq' — )• qq'Hj, mediated by the ex- 
change of a or a Z boson in the t-channel. The Higgs boson is radiated off this weak 
boson, while the photon can be emitted either from a fermion line or from a t-channel 
exchange W boson. The relevant charged-current (CC) Feynman diagrams can thus be 
grouped in two topologies, depending on how many gauge bosons couple to a fermion line. 
Representative graphs for each topology of a specific subprocess are depicted in Fig. |l[ 
Since, within the SM, the photon does not couple directly to Z bosons, only graphs corre- 
sponding to topology (b) contribute to neutral current (NC) production modes. 

For the calculation of the matrix elements we employ the methods applied already to a 



variety of WBF reactions pi, [2^, 26, 27, 2S], based on the helicity-amplitude techniques 



of Ref. [30, 31]. We decompose all Feynman diagrams into fermionic currents for each quark 
line and bosonic tensors parameterizing the gauge-boson interaction in the t-channel. For 
the graphs of topology (a), quark currents for both fermion lines and a bosonic tensor for 
the W~^W~ —7- Hj sub-amplitude are needed. Topology (b) furthermore requires a quark 
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Figure 1: Feynman diagrams contributing to the Born process us — >■ dcHj. Graphs analogous 
to (b), with the photon being emitted off the lower quark line, are not shown. 



current including photon emission and bosonic tensors for the sub-amplitude W~^W~ — )• H 
in CC processes and ZZ H in NC modes. We have developed three different imple- 
mentations of the tree-level matrix elements, which are supplemented by code adapted 



from (i) FeynArts |32] and Formcalc [33| and (ii,iii) HELAS |34|. In each implementation, 
building blocks entering in various diagrams are stored and evaluated once only. Diagrams 
for the related processes qq' — qq'H^, qq' — )• qq'H^, and qq' — )• qq'H^ are easily obtained 
by crossing. 

In addition to the topologies discussed above, annihilation processes such as qq — >• 
ZH^ with subsequent decay Z ^ qq and similar WH'y production modes occur. In 
sub-processes with identical quarks, interference contributions of t-channel with u-channel 
diagrams can arise. In the phase space regions where WBF processes can be observed 
experimentally, however, with widely separated quark jets of large invariant mass, these 
types of contributions are entirely negligible [^]. We therefore disregard them throughout. 
For simplicity, we refer to pp — )• Hjjj within the approximations discussed as "EW Hjjj" 
production. 

2.2 Calculation of the NLO-QCD corrections 

The calculation of real-emission corrections to EW H^jj production proceeds along the 
same lines as the LO computation. A gluon has to be attached to the fermions in all 
possible ways, yielding qq' — )• qq'gH^ contributions and related sub-processes with anti- 
quarks. Crossing the gluon to the initial state gives rise to (anti-)quark-gluon initiated 
modes with an additional fermion in the final state, such as gq' — )• qqq'H'j. The form 
of the bosonic tensors that already occurred at LO is not affected. Singularities in soft 
and collinear regions of phase space are handled in two different regularization schemes: 
in conventional dimensional regularization and in dimensional reduction, with d = 4 — 2e 
space-time dimensions. The cancellation of these divergences with respective poles in the 
virtual corrections is performed by introducing the counter-terms of the dipole subtraction 
formalism of Ref. IM]. Since the QCD structure of EW Hjjj production is identical to 
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the related case of pp — ?■ Hjj via WBF, the counter-terms are of the same form and can 



straightforwardly be adapted from Ref. |37], where also explicit expressions for the finite 
cohinear terms are given. 

The virtual corrections comprise the interference of one-loop diagrams with the Born 
amplitude. Due to color conservation, only selfenergy, vertex, and box corrections to either 
the upper or the lower quark line need to be considered. Contributions from graphs with 
a gluon being attached to both the upper and the lower quark lines vanish at order a^, 
within our approximations. The interference of the relevant diagrams containing one-loop 
corrections to either the upper or the lower quark line, My , with the Born amplitude, 
_B , is of the form 
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where Qi is related to the momentum transfer between the respective initial- and the 
final-state quarks carrying the momenta ki and k2 via Qf = — (fci — ^2)^, is the renor- 
malization scale, Cp = 4/3, and Cvirt is a constant, given by Cvirt = vr^/3 — 8 in conventional 
dimensional regularization and by Cvirt = 7r^/3 — 7 in dimensional reduction. The quantity 

■ (i) 

Aiy is a finite remainder. 

- — (*) 

In order to compute the M.v > we have split off the divergent pieces and expressed 
the remainder in terms of the finite parts of the Passarino-Veltman Bq, Cjk, and Djk 
functions which are evaluated numerically. To this end, we have prepared three 

diff'erent implementations, which are based on the Passarino-Veltman tensor reduction 
procedures of Ref. Ref. |3|, |4|, |4|], and Refs. [||, |4|, respectively. The tensor- 
reduction procedure may give rise to numerical instabilities for certain phase-space points, 
due to small Gram determinants emerging in the determination of the box-type corrections. 
We monitor these numerical instabilities carefully by requiring electroweak Ward-identities 
for all box-type contributions to be fulfilled with an accuracy of 10^'^ or better. We find that 
less than 0.1 permille of the generated events fails to meet this condition. The problematic 
contributions from these phase-space points to the finite parts of the box diagrams are 
negligible and therefore disregarded. 

The poles in Eq. ( |2.lD are canceled by respective singularities in the phase-space inte- 
grated counter-terms, which in the notation of Ref. are given by 
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2.3 Checks 



All building blocks entering the LO and NLO cross sections for EW H'jjj production 
have been tested extensively. As mentioned above, we have prepared three different im- 
plementations of the tree-level, real emission, and virtual matrix elements. We found 
perfect agreement between these implementations at the amplitude level. The tree-level 



- 4 - 



and real-emission contributions have also been compared to fully automatically generated 
amplitudes provided by MadGraph The matrix elements agree to about 12 digits for a 
representative set of phase-space points. 

Furthermore, we have verified the QCD gauge invariance of the real-emission contribu- 
tions, and the QED gauge invariance of the tree-level, virtual, and real-emission amplitudes 
by checking that they vanish upon replacing the polarization vector of the gluon and of 
the photon, respectively, with the corresponding momenta. 

In addition, we compared integrated cross sections at NLO-QCD accuracy within var- 
ious settings, including those of Sec. ^. All results agreed within the respective relative 
statistical errors, which are at the level of less than 10^^ for Monte-Carlo runs with sam- 
ples of about 50 million phase-space points. Integrated LO cross sections within inclusive 



selection cuts have also been compared to the corresponding results of MadEvent |46, 47 1. 
The cross sections agree within the accuracy of the two programs. 

A precise comparison of LO cross sections with those of Ref. |2^] was difficult, since 
not all EW parameters of the calculation have been listed in the publication. We have thus 
used our "default" settings of Sec. |^ for all parameters not explicitly listed in Ref. pOj . 
With this prescription, we could reproduce their cross sections with an accuracy of 3 — 5%, 
dependent on the chosen cuts. 

3. Numerical results 

For our numerical analysis, we assume a hadronic center-of-mass (cm.) energy of ^/S = 
14 TeV, unless stated otherwise. We use the CTEQ6M parton distribution functions with 
ctsiiTT-z) = 0.118 at NLO, and the CTEQ6L1 set at LO [^] as a default. We chose mz = 
91.188 GeV, mw = 80.398 GeV, and the measured value of Gp = 1.166 x lO'VCeV^ 
as electroweak input parameters, from which oqed and sin^ 6w are computed via LO 
electroweak relations. Throughout our calculation, fermion masses are set to zero, and 
contributions with external top-quarks are disregarded. Subprocesses with b quarks in the 
initial state are considered for the NC modes, while CC processes comprising the transition 
of a 6 to a t quark are not taken into account. For the Cabibbo-Kobayashi-Maskawa matrix, 
^KM) we have used a diagonal form, equal to the identity matrix, which is equivalent to 
employing the exact Vckm when the summation over final-state quark flavors is performed 
and quark masses are neglected. 



In order to reconstruct jets from the final-state partons, the algorithm [49, |50|, 51| 
as described in Ref. [^] is used, with resolution parameter D = 0.7. Jets are required to 
have 

pri>20GeV, Ivjl < 5 . (3.1) 

Here pTj denotes the transverse momentum of a jet, and yj the rapidity of the (massive) jet 
momentum which is reconstructed as the four-vector sum of massless partons of pseudo- 
rapidity \r]j\ < 5. At LO, there are exactly two massless final state partons, which are 
identified as tagging jets, provided they pass the kr algorithm and the cuts described 
above. At NLO, a third parton may be encountered which can either be recombined with 
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another parton or give rise to an additional jet. In this case, we choose to identify the two 
jets of highest transverse momentum as "tagging jets" . 

For the Higgs boson, we simulate a generic decay into two massless particles without 
specifying a particular channel. The decay particles, each one labeled d, can represent, for 
instance, bb final states. The respective branching ratio, BR{H — t- dd)^ is not included in 
the numerical results presented below. The photon is isolated in a theoretically well-defined 
way with the help of the criterion suggested in Ref . , which allows us to avoid introducing 
parton-to-photon fragmentation contributions. An event is considered as acceptable, if the 
hadronic energy deposited in a cone around the direction of the photon is limited by 

V PT^<^^^^^PT, (VAi2<Jo). (3.2) 

. . 1 — COS 00 

Here, the summation index i runs over all final-state partons found in a cone of size 
AR in the rapidity-azimuthal angle plane around the photon, pxi denotes the transverse 
momentum, and ARi^ the separation of parton i from the photon, while 5q stands for a 
fixed separation. 

In order to explore the impact of NLO-QCD corrections in different regimes, we have 
performed phenomenological analyses with two sets of selection cuts in addition to the 
common jet-defining criteria of Eq. ( |3.lD , to which we refer as "inclusive cuts" and "WBF 
cuts", respectively. 

For the "inclusive cuts" scenario, the Higgs decay particles are required to be separated 
from each other and from the other final-state particles by 

ARm > 0.4 , ARd^i > 0.4 , ARjd > 0.4 , ARj^ > 0.4 , (3.3) 

where each AR^k' (ARjk') denotes the separation of particle k (jet j) from particle k' in 
the rapidity-azimuthal angle plane. 

The cone-size parameter of Eq. (|3.2|) is set to 



6o = 0.7 . (3.4) 

The decay particles and the photon need to be located in the central-rapidity range of the 
detector, 

\yd\ < 2.5 , \y^\ < 2.5 , (3.5) 
and exhibit sufficiently large transverse momenta, 

PTd > 20 GeV , > 20 GeV . (3.6) 

The two tagging jets are required to fulfill the invariant mass criterion of 

Mjj^ > 100 GeV . (3.7) 

More stringent constraints are applied for the "WBF cuts" scenario. In addition to 



the cuts of Eqs. ( |3.1[ ), (3.5), and ( |3.6| ) we now require 

ARdd > 0.7 , ARd^/ > 0.7 , ARjd > 0.7 , ARj.^ > 0.7 , 6o = 0.7 . (3.8) 



-6- 



32 - 



28 



26 



r . Mo' = Q' (a) : 


.■'I 1 . . 1 . . ..| 1 i i 1 i i ii| . 
- (b) Aio' = niH'+ SpTj^ - 


LO fp= f 

NLO ?B= ?F= f 
NLO |p= f, f„= 1 
- NLO |r= f, fp= 1 

1 . . 1 ....1 1 . . 1 ....1 • 


■..1 1 . . 1 . . ..1 1 . . 1 . . ..1 ■ 



- 32 



28 



26 



0.1 0.2 0.5 1.0 2.0 5.0 10.0 0.1 0.2 0.5 1.0 2.0 5.0 10.0 



Figure 2: Scale dependence of the integrated cross section witliin the inclusive cuts of Eq. (|3.l| 



and Eqs. (3.3)-(B.7) at LO and NLO for two different choices of fig. Shown are curves for £,f = £, 
at LO (black dots), = = C at NLO (red solid), = = 1 at NLO (blue dash-dot), and 

^1 = 1 at NLO (green dashes). 



Backgrounds to WBF are significantly suppressed by imposing a large rapidity separation 
on the two tagging jets, 

^Vjj = IVji -yj2\ > 4 , (3.9) 

with the photon and the decay products of the Higgs boson being located in between the 
tagging jets, 

min(yji, yj2) < y-y, Vd < max(yji, ^^2) • (3.10) 
The tagging jets are furthermore required to reside in opposite detector hemispheres with 

2/jixyj2<0, (3.11) 

and exhibit a large invariant mass, 

M*J^> 600 GeV. (3.12) 

In order to estimate the dependence of our predictions on unphysical scales, we have 
computed the integrated cross section within the inclusive cuts of Eq. (3^) and Eqs. (3.3)- 
(3.7), fj'^'^*^, for two different choices of the factorization and renormalization scales, //p and 
/ipt, which are taken as multiples of the scale parameter /xq, 

/^f = Cf/^o, /^r = Cr^o- (3.13) 

Figure |2| (a) shows our results for = Qf, where for each fermion line //q is determined as 
the momentum transfer carried by the virtual weak boson emitted from it. In Fig. § (b), 
^cuts jg given for /Xq = mjj + ^ • , where for each event the sum runs over the transverse 
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^WBF[fb] for VS = 14 TeV 


LO 




CTEQ6 


MSTW 








f^l = Q\ 


^^l = m]J + Y.Plj 


0.5 


15.72 


14.56 


15.53 


14.30 


1.0 


14.65 


13.61 


14.40 


13.30 


2.0 


13.70 


12.76 


13.40 


12.42 




NLO 




CTEQ6 


MSTW 




/^g = Ql 




i^l = Ql 


^ll = m]J + Y.Plj 


0.5 


14.60 


14.84 


14.70 


14.93 


1.0 


14.79 


14.84 


14.91 


14.95 


2.0 


14.83 


14.75 


14.94 


14.85 



Table 1: Cross sections obtained for different values of the scale factor ^ = = ^i? within the 
"WBF cuts" scenario of Eq. (3.1) and Eqs. (5.S)-(3.12). The relative statistical errors of the quoted 
results are at the sub-permille level. 



momenta of all identified jets. Qualitatively, the results resemble the scale dependence of 
related WBF reactions [44|. In the range 1/2 < .^f = < 2, the LO cross section decreases 
by about 8% in each case, while the NLO prediction changes by about 2% for /Xg = 
and less than 1% for /Xq = rn^j + YIPtj- "^^^ ^ factor, defined as the ratio of the LO cross 
section to the respective NLO result, is close to one for the former choice (K = 1.02), while 
K = 1.07 for f^l = ml + ^p^^.. 

In order to quantify the impact of the parton distribution functions on the integrated 
cross section within the WBF-specific selection cuts of Eq. (3T) and Eqs. ( |3.8| )-( ^T^ ), 
^WBF^ in Table |l| we list the respective LO and NLO predictions as obtained with our 
default set, CTEQ6, and with the MSTW parton distributions of Ref. |Q for the two scale 
settings discussed above. Whenever the MSTW set is used, the corresponding expression 
for the strong coupling is employed. The differences between the LO cross sections for 
different parton distribution functions (but apart from that identical settings) are at the 
level of 2% and thus much smaller than those caused by the choice of the factorization scale. 
For instance, o"J^^(/i| = Qf) and <7Lcf^(//p = m'jj + YIPtj) differ by more than 7% for 
CTEQ6 and 8% for MSTW. In analogy to the "inclusive cuts" scenario discussed above, the 
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Figure 3: Rapidity separation [panel (a)] and invariant mass distribution of the two tagging jets 
[panel (b)] in EW Hjjj production at the LHC with ^/S = 14 TeV at LO (dashed black lines) and 



NLO (solid red lines), after the cuts of Eq. (3.1) and Eqs. (3.3)-(3.7) are applied 



scale dependence of a^^^ is mitigated by the inclusion of NLO-QCD corrections. For the 
NLO cross sections, the scale uncertainty is small and comparable in size to the uncertainty 
due to the parameterization of the parton distributions functions. In the following, we will 
use CTEQ6 parton distributions and set fi^ = Qf, unless stated otherwise. 

WBF-type reactions are characterized by widely separated hard jets in the far- forward 
and backward regions of the detector, being reflected by a large rapidity separation and 
invariant mass of the tagging jets. Figure |3| (a) illustrates the rapidity separation of the 
two tagging jets for the EW Hjjj cross section within the inclusive cuts of Eq. (|3.l| ) and 



Eqs. (U)-(U). Similar to the case of Hjj g^, |5|| and Hjjj [|§ production via WBF, 
the NLO-QCD corrections shift the peak of da/dAyjj to slightly larger values. Due to the 
possible presence of a third jet in the real-emission contributions, at NLO an enhancement 
of events with small values of Ayjj occurs. Such contributions can be efficiently removed 
by imposing the rapidity-separation criterion of Eq. (^^). The shape of the invariant 
mass distribution, depicted in Fig. ^ (b), is barely affected by the NLO-QCD corrections. 
Since da/dMj^^ peaks at rather large invariant masses and the distribution falls rather 
slowly at higher values of AfjJ^, the additional stringent cut of Eq. (3.12) is powerful in 
suppressing QCD backgrounds which exhibit invariant mass distributions with a much 
steeper slope than the WBF-type signal process. In the following we will therefore adhere 
to the WBF-specific cuts of Eqs. (|3.8|) ~( ^J2 ) in addition to the generic requirements of 
Eqs. (P), (ill), and (H). 

For this setting, the distribution of the hardest tagging jet is depicted in Fig. |^ (a). In 
order to assess the impact of the NLO-QCD corrections on the distribution of an observ- 
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Figure 4: Transverse momentum distribution of tlie hardest tagging jet in EW H^jj production 
at the LHC with ^/S = 14 TeV and \/S = 7 TeV, respectively, at LO (dashed black line) and 
NLO (solid red line) [panels (a) and (c)] and relative corrections according to Eq. ( 3.14 ) when 
the factorization and renormalization scales are varied in the range Qi/2 < /xr = fJ-F '^Qi 
[panels (b) and (d)]. 



able O, da/dO, together with the scale uncertainties of the LO and the NLO prediction, 
we consider the quantity d{0), defined as 



da{^F,^K)/dO 
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Figure 5: Transverse momentum distribution of the photon in EW H-fjj production at the LHC 
with ^/S = 14 TeV at LO (dashed black hne) and NLO (sohd red hne) [panel (a)] and relative 
corrections according to Eq. ( ^.14[ ) when the factorization and renormalization scales are varied in 
the range Qi/2 < fi^, = < 2Qi [panel (b)]. 



where dcr{S,F,S,ji)/dO denotes the LO or NLO expression, evaluated for arbitrary values of 
the scale parameters and ^r. The choice of /io is identical for da/dO and da^^^ /dO. In 
Fig. ^ (b), 6{p^j^) is shown for /Xq = and two different values of the scale parameters, 
^ = = = 1/2 and 2. The difference between the curves for the two values of ^ 
indicates the scale uncertainty of da/dp^j^ at LO (dashed black lines) and NLO (solid red 
lines), respectively. For low transverse momenta, the NLO-QCD corrections are positive 
and modify the LO results by more than 10%. With increasing pxj, the scale uncertainty 
of the LO prediction becomes large, amounting to about 19% for pTj = 200 GeV. The NLO 
prediction, on the other hand, is stable against scale variations over the entire transverse- 
momentum range considered, changing by less than 4% when ^ is varied from 1/2 to 2 
even for pxj = 200 GeV. Figure ^ (c) displays da / dp^r^^'^^^ for a hadronic cm. energy 
of ^/S = 7 TeV. While the size of the cross section obviously goes down with v^, the 
peak structure of the transverse momentum distribution is barely affected when the energy 
is decreased from 14 TeV to 7 TeV. The relative scale uncertainty of the LO prediction, 
illustrated by Fig. ^ (d), is significantly larger for a lower collision energy, however. 

In contrast to the hardest tagging jet, the photon exhibits a transverse momentum 
distribution rather insensitive to NLO-QCD effects. For our default choice, = A*f ~ Qi^ 
radiative corrections modify the LO result by less than 3% over the entire range of pT-y 
considered. The residual scale variation of the NLO-QCD prediction is very small. Our 
results for da/dpr-y and 6{pT'y) are shown in Fig. ^. 
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Figure 6: Distribution of tlie azimuthal angle separation of tlie two tagging jets in EW Hjjj 
production at the LHC with \/S = 14 TeV at LO (dashed black line) and NLO (solid red line) 



[panel (a)] and K factor according to Eq. (3.15) [panel (b) 



An observable particularly sensitive to the tensor structure of the Higgs coupling to 
weak bosons in WBF-type reactions is the azimuthal angle separation, A(pjj = \(j)ji — (j)j2\, 
of the two tagging jets ||7|, |]. Figure | illustrates da/dA(l)jj for "WBF cuts" with = Qf 
and ^/S = 14 TeV together with the phase-space dependent K factor, defined according to 



da^'^{lj,F,tJ-R)/dA(j)jj 



KiA<p,,) = - r^"^"^ iri" , (3.15) 



with the LO and the NLO distributions being evaluated for the same choice of scales. The 
shape of da/dAipjj is rather insensitive to NLO-QCD corrections. Should azimuthal angle 
correlations very different from this prediction be observed in experiment, they could thus 
hint at coupling structures not accounted for within the SM. 

The invariant mass distribution of the Higgs boson-plus-photon system, being recon- 
structed from the four-momenta of the photon and the decay products of the Higgs boson, 
is shown in Fig. |^ for our default settings with = Q? and ^/S = 14 TeV. The distribution 
vanishes for Mn-y < mn and peaks for mn = 120 GeV at around Mn^y ~ 165 GeV. For 
larger Higgs masses, da/dMn-y is shifted to correspondingly higher values. Reducing 
from 14 TeV to 7 TeV does not change the shape of the distribution, but increases the 
scale dependence of the LO prediction, in analogy to what has been observed above for the 
transverse momentum distribution of the tagging jets. 

4. Summary and conclusions 

In this work we have presented NLO-QCD corrections to Higgs production in association 
with a photon via weak-boson fusion at the LHC. We have developed a flexible parton-level 
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Figure 7: Invariant mass distribution of tlie Higgs boson-plus-photon system in EW Hjjj produc- 
tion at the LHC with ^/S = 14 TeV at LO (dashed black line) and NLO (solid red line) [panel (a)] 



and relative corrections according to Eq. (3.14) when the factorization and renormalization scales 
are varied in the range Qi/2 < /iR, = /ip < SQi [panel (b)]. 



Monte Carlo program which allows us to compute cross sections and kinematic distributions 
within experimentally relevant selection criteria, employing the photon-isolation procedure 



of Frixione |23] 



We analyzed EW H'jjj production within two different settings. First, we imposed 
only minimal selection cuts to obtain a well-defined final-state configuration. Second, 
additional cuts were applied, designed to enhance WBF-type contributions with respect 
to QCD background processes. We found that in each case the impact of NLO-QCD 
corrections on integrated cross sections is small. The actual size of the K factor depends 
not only on the selection cuts, but also on the choice of the factorization and renormalization 
scales in the computation. We studied the two options /^q = and /Ug = + YIPtj^ 
and found that slightly smaller NLO-QCD corrections are obtained for the former choice. 
The change of the NLO cross section when /uf and /iR are varied in the range ^o/2 < fip = 
A*R ^ 2^0 is comparable in the two cases. NLO-QCD corrections do not only affect the 
overall normalization of the integrated cross sections, but also the shape of some kinematic 
distributions. Relative corrections can be as large as 20% in some regions of phase space. 
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